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ABSTRACT 

The attenuation of very high energy 7-rays by pair production on the Galactic intersteUar radiation 
field has long been thought of as negligible. However, a new calculation of the interstellar radiation 
field consistent with multi-wavelength observations by DIRBE and FIRAS indicates that the energy 
density of the Galactic interstellar radiation field is higher, particularly in the Galactic center, than 
previously thought. We have made a calculation of the attenuation of very high energy 7-rays in 
the Galaxy using this new interstellar radiation field which takes into account its nonuniform spatial 
and angular distributions. We find that the maximum attenuation occurs around 100 TeV at the 
level of about 25% for sources located at the Galactic center, and is important for both Galactic and 
extragalactic sources. 

Subject headings: Galaxy: general — gamma-rays: observations — gamma-rays: theory — radiation 
mechanisms: general — radiative transfer 
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1. INTRODUCTION 

The attenuation of very high energy (VHE) 7-rays 
by pair production on the Galactic interstellar radia- 
tion field ( ISRF) has previouslv been considered to be 
negligible ijNikishovl IT961 lProtheroelll986|) . The main 
contribution is thought to come from pair production 
on the cosmic microwave background (CMB) where the 
effective threshold for attenuation is ^100 TeV and a 
maximum is reached at about 2000 TeV, currently ac- 
cessible only via air-shower experiments. The Galactic 
ISRF photons are more energetic so that the effective 
threshold is lower (^100 GeV) and the attenuation in- 
creases slowly to a maximum around 100 TeV. This cov- 
ers the energy range of present day Imaging Atmospheric 
Cherenkov Telescopes, such as the HESS instrument. A 
rough estimate of the attenuation of VHE 7-rays com- 
ing fr om the Galactic center (GC), which uses a new 
ISRF porter fc StroTigl.2nn5^ . but assumes an isotropic 
angular distri bution for the ISRF , shows that the effect 
is observable ijZhang et alJl2005j) . Given the essential 
anisotropy of the Galactic ISRF with most of the pho- 
tons going outwards from the inner Galaxy, the effect 
depends on the position of the source of VHE photons 
and its orientation relative to the observer in the Galactic 
plane. We calculate the attenuation of VHE 7-rays due 
to pair production with the Galactic photon field using 
the total ISRF over the entire Galaxy on a fine spatial 
grid which takes into account the nonuniform spatial and 
anisotropic angular distribution of background photons. 



^ Also Kavli Institute for Particle Astrophysics and Cosmology, 
Stanford University, Stanford, CA 94309 



2. INTERSTELLAR RADIATION FIELD 

The essential ingredients to calculate the Galactic 
ISRF are a model for the distribution of stars in the 
Galaxy, a model for the dust distribution and proper- 
ties, and a treatment of scattering, absorption, and sub- 
sequent re-emission of the stellar light by the dust. We 
briefly describe our ISRF calculation, wh ich is a further 
development of the work reported by iPorter &: Strong 
(200^; full details will be given in a forthcoming paper 
(Porter & Strong, in preparation). 

Our stellar model assumes a type cl assification based 
on th at used in the SKY model of iWainscoat et alJ 
l)1992() . It includes 87 stellar classes encompassing main 
sequence stars, AGB stars and exotics. For each stellar 
type there is a local star number density, scale height 
above the plane, fraction of local number density in each 
of several discrete spatial components, and spectrum in 
standard photometric filters. The stars are distributed 
in seven geometrical components: thin and thick disc, 
halo, bulge, bar, ring, and spiral arms. Spectra for nor- 
mal stars are t aken from the synthetic spectral library of 
iGirardi et al.l (|2.002) . Spectra for AGB stars and exotics 
are as given in the SKY model. 

We assume a dust model including graphite, poly- 
cyclic aromatic hydrocarbons (PAHs), and silicate. Dust 
grains in the model are spherical and the absorption 
and scattering efficiencies for graphite, PAHs, and sil- 
icate grains are taken from iLi fc Draina (|2001il . The 
grain model abundance and size d istribution are taken 
from IWeingartiier fc Draind l)200H) (their best fit Milky 
Way model), and a purely neutral interstellar medium 
is assumed. We consider only coherent scattering, 
and a Henyey-Greenstein a ngular distribution function 
ijHenvev fc Greensteinll94l(l is used in the scattering cal- 
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Fig. 1. — Interstellar radiation field energy density. Local in- 
terstellar radiation field (upper): heavy solid line, total radiation 
field including CMB; heavy dashed line, contribution by stars 
heavy dotted line, scattered light; heavy d ot-dashe d line, infra- 
red; thin solid line, local ISRF from lStrqng j^ Moskalenko fc Reime: 

§00). Data: squares , Apollo ^He nrv. Anderson & Fasti< 
10); triangles, DIRBE (l Arendt. et aTF iQQS): circles, FIRAS 
nkbcincr. Davis & Schlcgcl 1999). Interstellar total radiation 
field radial variation (lower): solid line, (R, z) = (0 kpc, kpc); 
dashed line, {R, z) = (4 kpc, kpc); dotted line, {R, z) = (12 kpc, 
kpc); dash-dotted line, (H, z) = (16 kpc, kpc). Heavy lines are for 
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our ISRF; thin lines for the ISRF of lStrong. Moskalenko 
J2000). 



culation. The stochastic heating of grains smaller than 
~ 0.1 /Ltm is treated u sing the "thermal continuous" ap- 
proach of IDraine fc Lil ( 2001); wc calculate the equilib- 
rium heating of larger dust grains by balancing absorp - 
tion with re-emission as described bv lLi fc Draind l)2001|) . 

Dust is assumed to follow the Galactic gas distribution. 
We use t he gas model for neutral a nd molecular hydrogen 
given by iMoskalenko et alJ l)2002|) . The radial variation 
in the G alactic metallicity g radient is taken to be 0.07 
dex/kpc ([Strong et al.ll20 (Mr and references therein). 

A cylindrical geometry is adopted for the radiation 
field calculation. Our calculations are simplified by 
assuming symmetry in azimuth and about the Galac- 
tic plane. The maximum radial extent is taken to be 
^max = 20 kpc. The maximum height above the Galac- 




FlG. 2. — Angles involved in the calculation: {R,z,a), galacto- 
centric coordinates of the source of VHE photons; {62,4>2), angles 
of the VHE photon; GC marks the Galactic center; and the ob- 
server's position is marked by "sun". 



tic plane is taken to be z^ax = 5 kpc. The Sun is located 
at Rs = 8.5 kpc from the GC. The Galaxy is divided into 
elements of equal volume and the total radiation field is 
calculated for each. 

The radiation field from stellar and scattered light is 
obtained using a modified form of the so -called partial 
intensity method l|Baes fc Deiongh3l200lD . Henceforth, 
wc will refer to the total stellar and scattered light as the 
"optical" radiation field. The infra-red radiation field is 
obtained by using the optical radiation field to calcu- 
late the dust emissivity for stochastic and equilibrium 
heating. The dust emission is integrated to obtain the 
infra-red radiation field throughout the Galaxy. 

Fig. n (upper panel) shows our calculated local ISRF 
including the CMB; also shown in the figure is the lo- 
cal IS RF calculated by iStrong. Moskalenko fc ReimeJ 
l)2000|) . The agreement of our computed ISRF with 
the data is generally good. We note that our new 
ISRF provides an improved fit to the observations 
around 10 /im c ompared t o the earlier model of 
IStrong. Mo skalenk o fc I^.eimeil (jToOG). This is particu- 
larly important for our optical depth calculation, since 
the majority of the attenuation is on these more numer- 
ous less-energetic photons of the ISRF. 

Fig. n (lower panel) shows the radial variation in the 
Galact ic plane of our ISRF (thick lines) to gether with 
that of lStrong. Moskalenko fc Reimeil l)2000() (thin lines) . 
Toward the inner Galaxy, our ISRF predicts a signifi- 
cantly higher energy density than previously described, 
particularly for the infra-red component. This arises be- 
cause of the coupling between the optical radiation field 
and the infra-red emission: in our calculation, the optical 
radiation field is used as direct input to the dust heating 
calculation, which was not done in previous work. The 
intense optical radiation field toward the inner Galaxy 
heats the dust to warmer temperatures, increasing the 
emission in the infra-red. Even though the optical emis- 
sion does not cause a direct enhancement to the attenu- 
ation at GeV to TeV energies - the energy of these pho- 
tons is typically a few eV, lowering the threshold, how- 
ever their number density is too low to provide significant 
absorption - it is essential to calculate this component 
of the ISRF to obtain the correct emission and angular 
distribution for the infra-red. Further discussion of the 
new ISRF is deferred to a forthcoming paper (Porter fc 
Strong, in preparation). 
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3. CALCULATIONS 

The optical depth for VHE 7-rays is given by the gen- 
eral formula: 



dN{e,ni,x) 



'11 



(ec)(l 



(1) 

COS 6'), 



where dN{e, fli,x)/ dedili is the differential number den- 
sity of background photons at the point x, e is the 
background photon energy, dfli — dcos9id4>i is a solid 
angle, cr^^ is the total cross section for the pair pro - 
duction process 77 e+e~ (I.Tau ch fc R.ohrlichlll980D . 
Ec = [■^eE{l — cos 6*)]^/^ is the center-of-momentum sys- 
tem energy of a photon, and 9 is the angle between the 
momenta of the two photons in the observer's frame. The 
integral over x should be taken along the path of the 7- 
rays from the source to the observer. 

The ISRF is cylindrically symmetric so that the pho- 
ton angular distribution depends on R and z only. In 
Fig. 121 we illustrate the galactocentric coordinate system 
{R,z,a). To calculate cos 9, the polar and azimuthal 
angles of the VHE photon, 92 and (l>2, are derived: 

= +RI-2RRs cos a, 
t>2 = -{Rs/p) sin a, 

-1/2 




(2) 



COS 



1 sin 02 cos((^i - 4>2), 



The 



where Rg is the galactocentric radius of the Sun. 
integration of eq. 1^ is done numerically. 

For the calculation of the optical depth in the CMB 
field we use the formula: 
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where kT is the CMB temperature, and rrieC^ is the elec- 
tron rest mass. 

4. RESULTS 

Fig. 131 shows the attenuation for selected positions 
(i?, z,a) as a function of incident 7-ray energy. For 
sources located at the GC the attenuation is ^12% at 
30 TeV and ~23% at 100 TeV. In Table d we give our 
optical depth results for selected values of R and z for 
a = 0°,90° and 180° at 30 and 100 TeV without con- 
tribution by the CMB. The attenuation is strongest for 
VHE 7-ray sources located toward the inner Galaxy and 
on its farside. 

To illustrate the effect of the anisotropic radiation field 
on the attenuation calculation, we show in Fig. 01 the 
ratio of the optical depths T^-y/ri^", where r-y^ is calcu- 
lated using the full angular distribution of the ISRF and 
t}^° is calculated assuming an isotropic distribution, for 
a source located at z = kpc (upper panel) and z — 5 
kpc (lower panel) emitting 100 TeV 7-rays as a function 
of position. The variation of t^^/t]^° over the Galaxy 
as seen from Earth is non-trivial. For sources located 
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Fig. 3. — Transmittance of VHE 7-rays as a function of 7-ray 
energy. Solid line: (i?, z,o) = (0 kpc, kpc, 0°) — L = 8.5 kpc; 
Dashed line: (-R, 2,0:) = (20 kpc, kpc, 90°) - L = 21.8 kpc; 
Dash-dotted line: {R,z,a) = (20 kpc, kpc, 180°) - L = 28.5 
kpc. Thick lines give the total transmittance curve including the 
ISRF and CMB. Left-most thin lines give the transmittance for the 
ISRF only; right-most thin lines for the CMB only. 



TABLE 1 





Optical depth 


Tjj AT 30 AND 100 TeV 






R, kpc 


z, kpc 


30 TeV 




100 TcV 




a = 0° 


90° 


180° 


a = 0° 


90° 


180° 








0.12 






0.25 






5 





0.01 


0.15 


0.22 


0.05 


0.32 


0.41 


10 





0.01 


0.08 


0.28 


0.02 


0.21 


0.54 


15 





0.02 


0.07 


0.31 


0.05 


0.18 


0.58 


20 





0.03 


0.07 


0.32 


0.06 


0.16 


0.60 





5 


0.03 






0.09 






5 


5 


0.01 


0.04 


0.06 


0.05 


0.10 


0.15 


10 


5 


0.02 


0.05 


0.09 


0.04 


0.12 


0.21 


15 


5 


0.03 


0.05 


0.12 


0.05 


0.13 


0.26 


20 


5 


0.03 


0.06 


0.14 


0.06 


0.13 


0.29 



in the Galactic plane between the GC and solar system 
the ratio is less than one since the majority of the ISRF 
photons are coming from the GC direction, cos — > 1 in 
eq. Q and the interactions are mainly following, leading 
to a lower pair production probability when the ISRF 
angular distribution is taken into account. The reverse 
situation, Tj^/t}^° > 1, occurs for 7-rays interacting with 
the ISRF in the outer Galaxy, where the majority of in- 
teractions are now head-on, cos 9 ^ — 1 in eq. For 
sources located beyond the GC, the ratio is '--^ 1. The in- 
terpretation of this case is straightforward: 7-rays emit- 
ted on the farside of the Galaxy toward the solar system 
have mainly head-on absorption interactions until they 
reach the GC, whereupon the interactions become mainly 
following. The isotropic case averages the angular distri- 
bution, and therefore leads to t^^/tI^° ~ 1. We note that 
the optical depth ratio increases monotonically with z, 
and is a result of the progressively more head-on nature 
of the 7-ray absorption interactions for sources located 
at larger distances from the plane. 



4 



Moskalenko et al. 




Fig. 4. — Effect of anisotropic ISRF on optical depth calculation at 100 TeV. Shown is the ratio of the optical depths T-y^/ri^" as a 
function of position {x,y) for a source located at 2 = kpc (left) and z = 5 kpc (right). The solar system is located at (x,y) = (8.5 kpc, 
kpc), exactly where the optical depth ratio rises sharply. 



Our results show that the attenuation of VHE 7-rays 
by the Galactic radiation field may be marginally ob- 
servable by the HESS instrument whic h has an effective 
sensitivity up to several tens of TeV ijAharonian et al.l 
I2QO5). The attenuation of VHE 7-rays from the sources 
on the Galaxy's farside will essentially steepen their spec- 
tra above ~10 TeV. In any case correction of source spec- 
tra for absorption is required. Interestingly, observations 
by future high energy experiments of the steepening of 



the spectra of Galactic sources in the GC region and be- 
yond may serve as a probe of the Galactic ISRF. 
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